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An  Iterative  Phase-Only  Nulling  Method 


1.  INTRODUCTION 

In  a  recent  report1  we  investigated  in  some  detail  the  problem  of  perturbing  a 
given  set  of  array  element  weights  to  impose  nulls  in  the  pattern  at  specified  loca¬ 
tions.  Two  types  of  perturbations  were  examined:  perturbations  of  both  the  ampli¬ 
tudes  and  phases  of  the  element  coefficients,  and  perturbations  of  the  phases  alone. 
When  both  the  amplitudes  and  the  phases  of  the  weights  are  allowed  to  vary,  the 
requirement  of  nulls  at  the  specified  locations  in  the  pattern  leads  to  a  system  of 
linear  equations  for  the  weight  perturbations  which  can  be  solved  exactly  subject  to  an 
additional  constraint  involving  minimization  of  the  perturbations.  If  only  the  phases  of 
the  element  coefficients  are  allowed  to  vary,  the  system  of  equations  for  the  imposed 
nulls  in  the  pattern  is  non-linear  and  cannot  be  solved  exactly.  Under  the  assump¬ 
tion  that  the  phase  perturbations  are  small,  however,  the  system  of  equations  for 
the  nulls  can  be  linearized  and  then  solved  in  the  same  way  as  the  equations  for  com¬ 
bined  amplitude  and  phase  perturbations.  Because  of  the  approximation  involved  in 
the  linearization,  the  resulting  phase  perturbations  do  not  give  perfect  nulls  as  do 
the  combined  amplitude  and  phase  perturbations.  The  depth  of  null  obtained  depends 
on  how  well  the  small  angle  approximation  used  to  linearize  the  equations  for  the 

(Received  for  publication  18  February  1982) 

1.  Shore,  R.  A.,  and  Steyskal,  H.  P.  (1982)  Nulling  in  Linear  Array  Patterns  With 
Minimization  of  Weight  Perturbations.  RADC-TR-82-32. 


nulls  is  actually  satisfied  by  the  solutions  to  the  linearized  equations.  It  was  found 
that,  even  for  a  low  sidelobe  array  in  which  the  linearized  phase-only  nulling  pro¬ 
cedure  gave  reasonably  deep  nulls  for  one  or  two  constrained  null  locations,  the 
effectiveness  of  cancellation  in  a  sector  of  the  pattern  achieved  with  phase-only 
nulling  deteriorated  as  the  size  of  the  sector  was  increased  by  adding  equispaced 
null  locations.  This  behavior  contrasted  sharply  with  that  found  when  combined 
amplitude  and  phase  perturbations  were  used.  The  combination  resulted  in  increas¬ 
ingly  effective  cancellation  throughout  a  sector  as  the  width  of  the  sector  was  in¬ 
creased  by  adding  more  nulls. 

The  failure  of  the  linearized  phase-only  nulling  procedure  to  give  a  satisfactory 
solution  to  the  nulling  problem  for  several  closely  spaced  nulls,  even  in  a  low  side- 
lobe  structure,  makes  it  of  interest  to  consider  alternate  procedures  for  solving 
the  phase-only  nulling  problem.  In  this  report,  1  describe  a  simple  scheme  of 
phase-only  nulling  in  which  the  linearization  procedure  is  used  iteratively.  This 
scheme  is  found  to  be  extremely  effective  in  low  sidelobe  structures  although  it  fails 
to  work  as  a  general  method.  The  method  is  described  in  the  next  section.  Some 
numerical  results  are  then  presented,  comparing  the  iterative  phase-only  nulling 
method  with  the  simple  (that  is,  non-iterative)  linearization  method  and  with  nulling 
using  combined  amplitude  and  phase  perturbations. 


2.  ANALYSIS 


Consider  a  linear  array  of  equispaced  isotropic  elements  (Figure  1).  The 
spacing  between  the  elements  is  d  and  the  phase  reference  center  is  assumed  to  be 
the  center  of  the  array.  Let  w  ,  n  =  1,2...,  N,  be  the  complex  weight  of  the  n**1 
array  element.  Then  the  array  field  pattern,  p(u)  is 


N 

p(u>  =  ^2 

n=  1 


w 

n 


u 


whe  re 


d  = 
n 


N  -  1 


-  (n  -  1),  n  =  1,  2, 


N 


and 


i 


u  =  kd  sin  Q 

' 

with 


6 


and  9  the  angle  measured  from  broadside  to  the  array.  Note  that  the  d^  are  odd 
symmetric  with  respect  to  the  phase  reference  center;  that  is, 

dn  =  -  dN  -  n  +  1  - 


Figure  1.  Geometry  of  Array 


The  general  problem  in  which  we  are  interested  is  as  follows.  Let 
an*  ns  1.2,  . . . ,  N,  be  a  given  taper,  assumed  to  be  symmetric  with  respect  to  the 
phase  reference  center,  of  the  amplitudes  of  the  element  excitations,  and  let  ug  be 
the  direction  for  the  peak  of  the  array  pattern.  Then  the  array  coefficients  are 


w  =  a  e 
on  n 


-id  u 
J  n  s 


n  =  1,2,  . . . ,  N  . 


We  wish  to  find  the  set  of  perturbations,  0^,  n=  1,  2,  . . . ,  N  of  the  phases  of  the 
element  weights  that  will  (a)  result  in  a  perturbed  pattern  with  nulls  at  a  set  of 
prescribed  locations,  u^*  k=  1,  2,  ...,  M,  and  (b)  be  "as  small  as  possible."  What 
does  it  mean  to  be  "as  small  as  possible?"  The  perturbed  coefficients  are 


w  = 
n 


a  e 
n 


■id  u 
J  n  s 


■j  d  u 
_  J  n  s 


■j  d  u 
e  n  3  <e 


j?n 


-  1)  . 


(1) 


The  first  term  on  the  RHS  of  Eq.  (1)  is  the  original  weight,  and  the  second  term  is 
the  total  perturbation  of  the  weight.  Following  Reference  1,  the  discussion  will 
consider  two  forms  of  minimizing  the  perturbations.  The  first  is  to  minimize 
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and  the  second  is  to  minimize 


l  (suJv^ 


The  first  form  corresponds  to  minimizing  the  perturbations  relative  to  the  original 
weights,  while  the  second  corresponds  to  minimizing  the  total  perturbations. 

From  Eq.  (1)  it  follows  that  we  can  represent  the  perturbed  array  pattern  as 
N 


-j  Uo  J  < 


n=  1 


a  e 
n 


(e  -  1)  e 


jd„u 


p(u)  =  p  (u)  +  ^  “  n  s  n  ,».n 


where  PQ(U'  is  the  original  array  pattern 


N 


-jdnU«  JdnU 


°  n  =  1  n 


Hence  the  equation  system  for  the  nulls  is 
N 


£  3n 


/idnusJ*n  nJdnuk 
e  (e  -  1)  e 


n  =  1 


po(uk),  k=  1.2.  ....  M  . 


(2) 


Separating  the  real  and  imaginary  parts  of  Eq.  (2)  we  obtain 
N 

L  an  cos[dn(us  -  uk)  -«nl  “  0 


n=  1 
N 


E  an  sinldn  <us  ‘  V  '^n1  =  ° 


n=  1 


k  =  1,2,  ....  M  . 


(3a) 


(3b) 


The  problem  can  then  be  stated  as  follows:  find  the  set  of  phase  perturbations, 
0n<  n-  1, 2 . N,  that  satisfy  the  equation  system  (3)  and  that  minimize  either 


n=  1 


.  2  Jn  . 
sin  (-j-) 


(4a) 


0  2 
1 11  n  sin  <‘T)f 


(4b) 
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This  problem  is  an  example  of  what  is  known  in  the  operations  research  litera¬ 
ture  as  a  nonlinear  programming  problem  (see  for  example  References  2  and  3). 

The  problem  in  general  has  no  exact  solution.  Furthermore,  it  is  not  at  all  obvious 
that  the  equation  system  for  the  nulls  has  any  solution  at  all  and  indeed  an  example 
can  easily  be  found  when  it  does  not  (consider  an  array  with  an  odd  number  of 
elements  for  which  the  amplitude  of  the  central  element  is  greater  than  the  sum  of 
all  the  other  amplitudes).  However,  we  will  not  attempt  here  to  answer  the  question 
of  under  what  circumstances  solutions  exist  to  Eqs.  (3).  Rather,  we  will  simply 
assume  that  the  given  amplitude  distribution  is  such  that  the  equation  system  for 
the  pn  is  underdetermined  if  2M  <  N,  that  it  has  an  infinity  of  possible  solutions, 
and  that  a  minimization  criterion  is  used  to  select  a  particular  one  of  these  solutions. 
Since  an  explicit  solution  to  the  problem  cannot  be  found,  an  efficient  numerical 
method  for  finding  the  solution  is  desired. 

The  scheme  we  propose  here  is  based  on  the  assumption  that  the  sidelobes  of 
the  original  pattern  are  sufficiently  low  that  only  small  phase  perturbations  are 
required  to  place  nulls  at  the  specified  locations.  Under  this  assumption  we  can 
expand  Eqs.  (3a,  3b)  and  neglect  all  terms  quadratic  and  higher  in  <i>n  thus  obtaining 


N 


an°l,nsinfdn(us  '  V  =  '  W 

n=  1 

N 

S  an°l,nCOS  l‘n  <Us  ~Uk)1=  ° 
n  =  l 


k  =  1,2,  ....  M  . 


(5a) 


(5b) 


Under  the  same  assumption  of  small  0  the  functions  to  be  minimized  can  be  approxi- 

n 

mated  by 
N 


n=  1 


(6a) 


and 


N 


(6b) 


We  have  written  o  instead  of  on  here,  not  only  to  indicate  that  the  solution  to 

Eqs.  (5)  is  an  approximation  to  the  solution  to  Eqs.  (3)  we  are  seeking,  but  also, 

as  will  be  seen  shortly,  because  the  O.  are  to  be  the  first  stage  in  an  iterative 

l » n 

process  leading  to  the  desired  solution. 

2.  Himmelblau,  D.  M.  (1972)  Applied  Nonlinear  Programming,  McGraw-Hill,  NewYork. 

3.  Schittkowski,  K.  (1980)  Nonlinear  Programming  Codes,  Springer-Verlag,  NewYork. 


The  solution  to  Eq.  (5a)  that  minimizes  Eq.  (6a)  or  Eq.  (6b)  was  found  in  our 
earlier  report1  and  the  result  is  merely  stated  here:  let 

n  ; 

2  ♦?.  n  =  min 


C  = 

n 


a  , 
n 


n  =  1 
N 


!•  E  <an^.n)2  =  min' 


n  =  1 


(7) 


n  N 


(8) 


n  =  1 


and 


1  =  -  lP0(Ul).  P0(u2).  P0(um)J  ; 


then 


2 1  ^ 

l.n  “  bm  Sin  [dn  (Us  -  Um)] 


(9) 


where  the  vector  of  coefficients  bm  is  obtained  by  inverting  the  matrix  equation 
A  b  =  £ 


with  the  elements  of  the  matrix  A  given  by 
N 


W=  ^2  2t  sin  [ci  (u  -  uk>]  sin  [d  (u  -  u  )]  k,  m=l,2 . 

km  n=  1  n  n  s  k  n  s  m 


M 


Equation  (5b)  is  automatically  satisfied  by  this  solution  because  of  the  antisymmetry 

of  the  0.  with  respect  to  the  phase  reference  center.  The  weights  corresponding 
i»  n 

to  the  phase  perturbations  <p  ^  R  are 


-jdnus  j0l  n 
w,  =  a  e  nSgi.n 

1,  n  n 


and  the  array  pattern  is  then 

V1  J  u 

?!<">  =  n?j  WJ,ne 
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Now,  if  the  phase  perturbations  <p ^  n  are  fairly  good  approximations  to  the 
desired  perturbations  0n>  then  the  values  of  p^(u)  at  the  specified  null  locations 
u  =  u^,  k  =  1,  2,  ....  M,  will  be  smaller  in  magnitude  than  the  original  values 
pQ(u^).  Hence  it  is  reasonable  to  suppose  that  the  coefficients  Wj  can  be  per¬ 
turbed  in  turn  with  a  set  of  still  smaller  phase  perturbations  so  as  to  place  nulls  in 

the  pattern  at  the  specified  locations.  This  idea  leads  us  to  the  following  iterative 

4  5 

scheme,  similar  to  that  used  by  Cheng  and  Raymond,  and  by  Cheng.  Let  the 
weights  at  the  end  of  the  i1^1  iteration,  i  -  1,  be 


i,  n 


a 


n 


-J 

e 


d 


n 


u 

s 


n 


n  =  1,2,  ....  N 


and  the  corresponding  array  pattern 

N  ■  ^ 

EJ  d  u 

w.  e 

n-  1  1*n 


For  i=0,  define 
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o. 


n 


0 


w 


o,  n 


a 


n 


e 


n 


u 

s 


JLet  <J>.+  1  ^  -  <J>.  n  be  the  phase  perturbation  introduced  at  the  (i+  l)**1  iteration, 
i  -  0.  Then  similarly  to  Eqs.  (5a  and  5b),  we  have 


N 


2 L*  a  (0.  ,  .  -  0.  )  sin  f  d 

"j  n  ri+  1,  n  Yi,n  1  n 


(us  -V  ~  *i.n>  =  "  PiV 


N 


2~i  a  (0.  ,  -  0.  )  cos  [d  (u  -  u.  )  -  0.  ]  =  0 

,  n  ri+  l,n  yi,  n  1  n  s  k  ri,  n' 
n  =  i 


k  =  1,2, 


M 


or 


4.  Cheng,  D.K. ,  and  Raymond,  Jr.,  P.  D.  (197 1)  Optimization  of  array  directivity 

by  phase  adjustments.  Electron.  Letts.  ^552-554. 

5.  Cheng,  S.  D.  (197  1)  Optimization  techniques  for  antenna  arrays,  Proc.  IEEE. 

90:1664-1674. 

vw 
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Za  0.  ,  sin[d  (u  -  u,  )  -  o.  1  =  -  p.(u.  ) 

j  nri+l,n  n  s  k  1.  nJ  k 

N 

+  X  a  o  sin  [d  (u  -  u,  )  -  q. 

“j  n  l,  n  1  n  s  k  i,  n 


(10a) 


£  <a  0.  )5 

n=i  n  l+1*n 


Eqs.  (10a  and  10b)  are  of  the  same  form  as  Eqs.  (5a  and  5b)  except  for  the  RHS  of 
Eq.  (10b).  But  because  of  the  antisymmetry  of  the  Oj  l  noted  above,  the  KHS  of 
Eq.  (10b)  is  zero  when  i=  1,  and  hence  by  induction  the  same  will  be  true  for  succeed¬ 
ing  iterations  as  well.  Thus,  similarly  to  Eq.  (9),  the  0.+  j  ;  are  given  by 

2 1 

0.  ,  =  — -  X  b.  sin  [d  (u  -  u  )  -  o.  ]  .  (11) 

l  +  1 ,  n  a  i+l.m  1  n  s  m  i,  nJ 

n  m  =  1 


The  vector  of  coefficients  b.  +  j  m  is  obtained  by  inverting  the  matrix  equation 
A.  .  b  =  y.  . 

i+  1  —i  +  1  —i+l 


where 


(12) 


Ei+l  =  'j  W"  2,  an^i.nsinIdn<us  ‘V-^nl . P^V 

N  T 

-  / *,  a  0.  sinfd  (u  -  U-.)  -  u.  ]  J 
nri,n  1  n  s  M  i,  n'  ) 


(13) 


?! 


Eao.  ,  cos  d  u  -  u.  )  -  o.  1=  X,  a  0  cos (d  (u  -  u.  )  —  ©.  ] 

,  r.  i+l,n  1  n  s  k  l.n1  *-*,  n  i,n  'ns  k  i,  nJ 
n=  1  n=l 

(10b) 


The  functions  to  be  minimized  are 
X 


E  of. , 


n=  1 


i  +  l.n 
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and  the  elements  of  the  matrix  A^+  ^  are  given  by 
N 

Ta.  ."I  =  E  2  t  sin[d  (u  -  u,  )  -  p.  ]  sin  d 

i  +  1  I,  n  n  s  k  ri,n'  n 

L  Jk,  m  n=  1 


(uo  -  u  ) 
s  m 


l.  n 


I  .  04) 


Testing  for  convergence  and  termination  of  the  iterative  procedure  can  be  based 
either  on  the  phase  perturbations  themselves  (for  example,  termination  when 

m*x  {  K+En'^.n'l  <4  orwhen  E  l^i+i.n-^i.n1  < £  >•  or  on  the 

values  of  the  beam  coefficients.  That  convergence  of  the  phase  perturbations  implies 
convergence  of  the  beam  coefficients  follows  easily  from  Eq.  (12)  in  conjunction  with 
Eqs.  (13)  and  (14),  since 

b.  =  A.  *  .  v. 

—l+l  i+  1  *-i+  1 


and  A.  +  ^  and  y^+  ^  approach  limiting  values  if  the  (k  n  converge  to  limiting  values 

V 

It  is  also  worth  noting  that  if  the  limiting  values  of  the  phase  perturbations  are 

small,  the  resulting  cancellation  pattern  can  be  interpreted  as  the  sum  of  pairs  of 

slightly  distorted  beams,  one  pair  for  each  imposed  null,  in  a  way  that  is  completely 

analogous  to  the  interpretation  of  the  non-iterative  phase-only  cancellation  pattern 

as  the  sum  of  pairs  of  beams  (see  Reference  1).  For  from  Eq.  (11),  we  see  that  if 

the  0-  converge  to  limiting  values  O  ,  then 
i*  n  n 


2 1 


M 


=  — —  y.  b  sin  [d  (u  -  u  )  -  o  ]  . 
n  a  '  m  1  n  s  m  n 


(15) 


n  m  -  1 


Also,  from  Eq.  (1),  if  0n  is  small  then  the  change  in  the  n  weight  can  be  approxi¬ 
mated  by 

-j  d  u  jo 

a  n  s  .  J  n  , . 

A  w  =  a  e  (e  -  1 ) 


da  0  e 
J  n  n 


-id  u 
J  n  s 


(16) 


so  that,  substituting  Eq.  (15)  in  Eq.  (16), 

E  'jdnus  (  ild  <u  _u  -jfdn(u=  -u  )  '  Pnl  > 

be  n  b  <  n  s  m  n  nsm  n  l 

_ ,  m  J  e  ' 


ms  1 


{  -j(d  u  +  p  )  —  j[ d  (2  u  -  u  )  -  p  1  > 
t  £  b  e  n  m  9n  -  e  n  3  m  n  ( 
n  m =  1  m  (  ' 


13 


and  hence  the  cancellation  pattern  is  given  approximately  by 


M 


N 


Ap(u)  *51  bm  53  tn 


m=  1  m  n=  1 


/dn(u  '  Um}  ‘  ^n)].  ejIdn(u 


|2us  -um}  >+  0n] 


If  the  phase  perturbations,  0n»  were  zero,  the  cancellation  pattern  would  thus  be 
represented  as  the  sum  of  M  pairs  of  beams,  one  member  of  each  pair  directed  at 
an  imposed  null  location,  and  the  other  member,  of  opposite  sign,  directed  at  the 
symmetric  location  with  respect  to  the  axis  of  the  mainlobe  of  the  original  pattern. 
Indeed,  this  is  the  approximate  representation  of  the  cancellation  pattern  for  the 
noniterative  phase-only  nulling  scheme  (see  Reference  1).  If  the  phase  perturbations 
are  non-zero  but  small,  it  is  to  be  expected  that  the  phase  perturbations  introduce 
a  small  amount  of  distortion  in  the  shape  of  the  component  beams  but  that  the  beams 
representation  of  the  cancellation  pattern  still  holds.  Hence,  for  small  perturba¬ 
tions,  the  limiting  values  of  the  beam  coefficients  obtained  with  the  iterative  nulling 
method  will  not  differ  greatly  from  the  values  of  the  coefficients  obtained  with  the 
noniterative  nulling  scheme.  As  with  the  non-iterative  method,  the  nulls  produced 
at  the  prescribed  locations  will  be  accompanied  by  a  comparable  raising  of  the 
pattern  at  the  symmetric  locations  with  respect  to  the  mainlobe  axis. 

If  the  iterative  procedure  converges,  the  limiting  values  of  the  0i  n  should 

satisfy  the  system  of  Eqs.  (3)  for  the  pattern  nulls.  If  the  limiting  values  are  small 

they  should  also  be  reasonable  approximations  to  the  values  of  the  <f>n  that  not  only 

satisfy  Eqs.  (3)  but  which  also  minimize  either  expression  (4a)  or  (4b).  Further 

work  is  needed,  however,  to  obtain  a  quantitative  estimate  on  how  close  the  limiting 

values  of  the  rt.  will  be  to  the  ideal  solution.  A  direct  comparison  of  the  values 
l ,  n 

obtained  by  a  more  powerful  non-linear  programming  method  with  the  values  ob¬ 
tained  here,  would,  of  course,  serve  this  purpose.  The  solution  obtained  by  the 
iterative  scheme  proposed  here  has  value  not  only  because  it  provides  a  useful 
approximation  to  the  desired  answer,  but  also  because  it  can  serve  as  a  starting 
point  for  further  refinement  by  a  more  powerful  non-linear  optimization  method. 

Such  methods  often  must  be  started  at  a  feasible  solution— that  is,  a  solution  to  the 
constraint  equations. 


*On  first  glance  it  might  appear  that  the  solution  to  Eqs.  (3)  that  minimizes,  say 

53  also  minimizes  (4a).  For  if  we  let  {  0n|  be  the  set  of  phase  perturbations  that 
satisfy  Eqs.  (3)  and  minimize  (4a),  and  }0nj  be  the  set  of  phase  perturbations 

that  satisfy  Eqs.  (3)  and  minimize  ^0^  then  the  j^[  and  J®n}  must  satisfy  the 
following  set  of  inequalities,  assuming  0<  |0  l-it/2,  0<  \\p  ^\  <  it  /  2: 
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3.  NUMERICAL  RESULTS 

A  computer  program  was  written  to  implement  the  iterative  nulling  scheme 
described  in  the  previous  section,  and  run  on  a  CDC  6600  computer.  All  computa¬ 
tions  were  performed  for  an  array  with  41  elements  with  half  wavelength  spacing. 
Unless  otherwise  stated,  the  original  or  unperturbed  pattern  corresponds  to  a 
40  dB  Chebyshev  taper  of  the  element  excitations.  Figure  2  shows  results  obtained 
when  the  iterative  program  was  run  to  impose  nulls  at  one  (15.  23°),  two  (15.  23°, 

15.  78°),  three  (15.  23°,  15.78°,  16.33°),  and  four  (15.  23°,  15.78°,  16.33°,  16.88°) 
locations.  We  have  plotted  the  maximum  of  the  values  of  the  perturbed  pattern  power 
at  the  locations  at  which  nulls  were  desired  (that  is,  the  least  deep  null)  vs  number 

of  iterations.  Figure  2a  gives  the  results  obtained  when  53  (a^  <Pn)2  was  minimized. 


E  ..»24>s  4  £  s1„24,<  E  f2?  E  *2. 

n=  1  2  n=  1  *  n* 1  n  n=  1  n 


If  in  general 
N 


E  *1’-  E  e 

n=  1  n  n  =  1  n  —  n-  1  2  n- 1 


sin2  <!», 


<t> 


as  it  does  when  N=  1,  it  would  follow  that  the  first  "<"  in  Eq.  <t>  must  be  strictly 
=  .  Unfortunately  this  is  not  true,  as  can  be  seen  by  letting  j  =  =0,n>  1 

<J»n  =  n  =  1,2,  ....  N.  We  then  have,  for  0  <  e  <1,  ^  ^/2  >  53  02  and 

sin2  (-j)  =  sin2  (1/2)  =  0.  230 


but 


N 


sin2  (—4==. )  >  N 


n-  1 


2Vn' 


(— 1_  )2  -  i 
2  y??  3  2yN 


>  0.  240 


for  N  >  1. 
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and  figure  2b  gives  the  results  corresponding  to  minimizing  The  starting 

point  of  all  curves  is  the  maximum  value  of  the  power  of  the  original  pattern  at  the 
null  locations.  The  value  for  one  iteration  is  the  least  null  depth  achieved  with  the 
simple  (that  is,  non-iterative)  linearized  nulling  procedure  that  coincides  with  the 
first  step  of  the  iterative  procedure.  We  see  that  in  all  cases  the  iterative  procedure 
gives  null  depths  (<  -250  dB)  that  greatly  exceed  the  null  depths  obtained  with  the 
non-iterative  nulling  procedure.  As  a  number  of  null  locations  increases,  however, 
more  iterations  are  required  to  achieve  the  same  depth  of  null.  There  is  also  a 
tendency  (as  is  shown,  for  example,  in  the  four-null  curve  of  Figure  2b)  for  the 
maximum  depth  of  null  achievable  by  the  iterative  procedure  to  slowly  decrease  as 
the  number  of  null  locations  increases. 


NUMBER  OF  ITERATIONS 


Figure  2.  Maximum  of  Perturbed  Pattern  Power  at  Imposed 
Null  Locations  for  1,  2,  3  and  4  Nulls  as  a  Function  of 
Number  of  Iterations.  40  dB  Chebyshev  amplitude  taper. 
Figure  2a:  s(a  0)2  =  min;  Figure  2b:  £0“  =  min 


Figure  3  shows  corresponding  results  obtained  starting  with  a  20  dB  Chebyshev 
taper  of  the  element  excitations,  when  the  iterative  program  was  run  to  impose  nulls 
at  one  (14.  70°).  two  (14.  70°,  15.28°),  three  (14.  70°,  1528°,  15.86°),  and  four 
(14.70°,  15.28°,  15.86°,  and  16.  44°)  locations.  We  see  that  for  both  minimization 
criteria  the  iterative  procedure  was  able  to  produce  deep  nulls  for  the  one  and  two 
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null  cases,  although  considerably  more  iterations  were  required  compared  with  the 
40  dB  calculations  of  Figure  2.  (Note  the  curious  feature  of  Figure  3b  in  which 
fewer  iterations  were  required  for  the  two-null  case  than  for  the  one-null  case.  ) 

For  three  nulls,  the  iterative  procedure  was  ineffective  for  the  case  of  minimizing 

/.  (a  <J>  )“  but  worked  well  for  the  case  mmizing  o^.  For  four  nulls,  the 

iterative  procedure  was  unable  to  achieve  any  nulling  for  either  minimization  cri¬ 
terion.  These  results  point  out  the  limitations  of  the  iterative  method  when  the 
sidelobes  of  the  original  pattern  are  relatively  high.  The  phase  perturbations  re¬ 
quired  to  produce  nulls  in  the  pattern  are  then  larger  in  magnitude  than  they  are  for 
low  sidelobe  patterns,  especially  when  multiple  closely  spaced  nulls  are  prescribed, 
with  the  result  that  the  iterative  procedure  (consisting  of  repeated  linearizations 
based  on  the  assumption  of  small  phase  perturbations)  fails  to  work. 


NUMBER  Of  ITERATIONS 


Figure  3.  Maximum  of  Perturbed  Pattern  Power  at 
Imposed  Null  Locations  for  1,  2,  3  and  4  Nulls  as  a 
Function  of  Number  of  Iterations.  20  dB  Chebyshev 
amplitude  taper.  Figure  3a:  £( a  <t>  )2  -  min; 
Figure  3b:  x<tr  -  min  n 


As  an  application  of  the  iterative  nulling  method  calculations  were  performed 
to  study  the  variation  of  cancellation  effectiveness  within  a  sector  of  the  pattern  as 


17 


the  width  of  the  pattern  is  increased  by  adding  equispaced  nulls.  As  in  our  earlier 
report1  the  power  cancellation  ratio  in  the  sector  A 8  =  8^-8  -  0  j  defined 


max 
8  £  A0 


max 
8  £  A  8 


tP«?>]2 

[Po<0)]2 


where  p  (8  )  is  the  original  pattern  and  p(8  )  is  the  perturbed  pattern.  In  Table  1  we 
have  tabulated  the  cancellation  ratios  for  various  sector  widths  and  different  nulling 
methods.  The  entries  in  the  i  n  row  of  the  table  are  the  cancellation  ratios  ob¬ 
tained  by  the  methods  listed  across  the  top  of  the  table  for  the  sector 

15.  23°  s  $  <  {8  „).  with  nulls  placed  at  15.  23°  and  at  the  locations  8  =  (8  „).,  j  -  i. 

^1  £  ] 

For  example,  the  entry  of  -68.  0  dB  in  the  fourth  row  for  iterative  phase-only  nulling 

V*  2 

minimizing  (a^  0^)  means  that  when  nulls  were  imposed  by  this  method  at  the 
locations  15.23°,  15.78°,  16.33°,  16.88°,  and  17.44°,  the  cancellation  ratio  in  the 
sector  15.23°  5  8  -  17.44°  was  -68.0  dB.  The  sequence  of  null  locations 
15.23°,  ...  18.  00°  is  equispaced  with  respect  to  sin0  .  For  comparison  purposes 
we  have  included  in  Table  1  not  only  the  cancellation  ratios  obtained  with  the  two 
iterative  phase-only  nulling  methods  but  also,  taken  from  Reference  1,  the  cancella¬ 
tion  ratios  obtained  with  the  non iterative  phase~only  methods  and  the  exact  combined 
amplitude  and  phase  perturbation  methods. 

The  outstanding  feature  of  Table  1  is  the  remarkable  improvement  in  cancella¬ 
tion  effectiveness  resulting  from  the  use  of  the  iterative  phase-only  methods  as  com¬ 
pared  with  the  non-iterative  phase  -only  methods.  For  six  nulls,  cancellation  ratios 
of  -82  dB  and  -78  dB  were  obtained  using  the  iterative  methods  as  compared  with 
a  -1  dB  for  the  non-iterative  phase-only  methods.  This  improvement  is,  of  course, 
attributable  to  the  fact  that  the  iterative  methods  are  able  to  produce  far  deeper  nulls 
at  the  prescribed  locations  than  are  the  non-iterative  methods,  with  the  result  that 
the  entire  pattern  within  the  sector  is  pulled  down  much  further.  Indeed,  the  iterative 
phase-only  methods  give  sector  cancellation  comparable  to  the  exact  combined  ampli¬ 
tude  and  phase  perturbation  methods  for  all  values  of  02  through  17.44°.  For 
8 ^  3  18-  00°,  the  iterative  phase-only  methods  give  slightly  less  effective  sector  can¬ 
cellation  than  do  the  combined  amplitude  and  phase  perturbation  methods.  This  is 
probably  attributable  to  the  fact  that,  as  noted  above,  the  depth  of  null  achievable  by 
the  iterative  methods  slowly  degrades  as  the  number  of  imposed  nulls  increases. 

For  =  18.00®,  that  is,  six  imposed  nulls,  the  shallowest  of  the  six  nulls  produced 

by  the  iterative  phase-only  methods  was  -238  dB  for  minimized  S  and 

-221  dB  for  Z^2. 
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Even  though  the  iterative  phase-only  methods  studied  in  this  report  are  not  as 
effective  in  certain  multiple  null  situations  and  in  high  sidelobe  applications  as  are 
combined  amplitude  and  phase  perturbation  methods,  nevertheless  it  is  important 
to  recognize  that  these  limitations  are  not  limitations  so  much  of  phase-only  nulling 
per  se  as  they  are  of  the  particular  iterative  scheme  proposed  here.  Based  on  the 
results  obtained  and  discussed  in  this  report,  it  would  appear  that  in  many  situations 
phase-only  methods  are  in  principle  as  effective  in  nulling  and  sector  cancellation 
as  are  combined  amplitude  and  phase  perturbation  methods.  Further  work  is  needed, 
however,  to  refine  the  numerical  procedures  involved  in  phase-only  nulling.  The 
principal  drawback  of  phase-only  nulling  is  that  the  nulls  and  sector  cancellation 
achieved  on  one  side  of  the  mainlobe  of  the  original  pattern  are  accompanied  by  a 
comparable  raising  of  the  pattern  at  the  locations  and  sectors  symmetrically  placed 
with  respect  to  the  mainlobe.  This,  rather  than  limited  nulling  effectiveness, 
appears  to  be  the  main  price  one  pays  for  using  phase-only  nulling. 


4.  CONCLUSIONS 

In  this  report,  an  iterative  numerical  method  is  described  for  calculating  the 
minimum  phase-only  perturbations  of  the  element  excitations  of  a  linear  array 
required  to  impose  nulls  in  the  pattern  at  prescribed  locations.  The  method  is  based 
on  repeated  linearizations  of  the  equations  for  the  imposed  nulls.  The  iterative 
scheme  is  found  to  be  extremely  effective  in  low  sidelobe  applications  when  the  re¬ 
quired  phase  perturbations  are  small,  but  fails  to  work  as  a  general  method  when 
the  required  phase  perturbations  are  large. 
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